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Abstract. Let p be a prime integer and Hg be a collection of complex positive 
definite symmetric g X g matrices r. Denote by pr the multiplication of t by 
p. In this note we describe an explicit process to obtain algebraic identities 
between theta functions with integral characteristics evaluated at r and pr. 
For g = I this produces modular equations between X{t),X{pt) where A(t) is 
the invariant associated with elliptic curve generated by t, described by the 
equation: = x{x — 1) (2; — A(ti)) . Consequently, if g > 1 The algebraic 
identities we obtain might serve as a higher dimensional generalization for the 
one dimensional modular equations. 
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inspiration. 

1. Introduction 

Let Ti be a complex number such that Im(Ti) > and is the lattice generated 
by {l,Ti} in C Let Ci = C/Z^, be the corresponding analytic elliptic curve. The 
algebraic equation of this curve is: 

(1) = x{x - 1) (x - X{n)) . 

A(ti) is the invariant corresponding to ri in this equation. 

Definition 1.1. Let p be a prime number. A modular equation of order p for A, 
is an algebraic equation between A(ri) and A(pTi). 

These equations appeared naturally in the theory of elliptic integrals because, 
such an equation describes algebraically the analytical multiplication by p on the 
lattice Zt-j^ . In equivalent more contemporary terms this equation describes A in- 
variant of curves C'2 and (j) : C'l '—^ C2 is an isogeny (finite homomorphism) of 
order p. Equations of this type have an important role in Galois theory of complex 
multiplications since, if Ci is a curve with complex multiplication then A(pTi) gen- 
erates interesting field extension of (5°''(ri). Another application of one dimensional 
modular equations is algorithms for rapid calculation of tt, [Bo . 

In recent years there are applications of modular equations to point counting 
algorithms of elliptic curve above finite fields Fpi . Modular equations of order p are 
used to compute explicit canonical lifting of elliptic curves above finite fields Fpi , 
of characteristics p to the corresponding p - adic field above it. Using the lifting we 
calculate the trace of the Frobenius operator on the p- adic field. Applying fixed 
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point formulas, we find the number of points of elliptic curves above finite fields 
Fpi. [Ma| explains the general framework for this type of algorithms. 

For p = 2 the modular equation of level 2 is the arithmetic geometric mean 
(AGM) in disguise. If ao = a, 6o = & are real numbers we define the AGM iteration 
as: 



(2) 



(3) 



bn = V Clnbn 



This iteration has a strong link with the modular equations of order 2. These se- 
quences converge to a common limit denoted by AGM{a,h) see |Bo) . Mestre |Me] 
used the AGM iteration to suggest an algorithm for point counting defined over 
F2i. |Me| Mestre uses a higher dimensional analogue of the AGM and produces 
an algorithm to count the number of points of hyper-elliptic curves above fields 
of characteristics 2. In a recent work Lubicz Carls and Kohel jCKLj generalized 
Mestre's method further in a different direction. They construct curves that ac- 
cording to them have good cryptographic properties using an iteration of order 3. 
In one dimensional case their iteration seems to be closely related to modular equa- 
tions of order 3. These results motivate the question whether there exists a theory 
of modular equations for r an element of Hg. To suggest a possible generalization 
recall that A(ti) is a quotient of theta functions that is: 



A(ri) = 



94 


■ ■ 

1 


(0,ri) 


64 


■ ■ 




(0,n) 



These are analytic functions that we define in the first section of the paper. 
A(pri) is the same quotient evaluated at pTi. Hence modular equations become 
identities between theta functions evaluated at ti and those evaluated at pri . While 
an elementary dimension argument shows that analogues of A(ri) do not exist for 
general r G Hg theta functions do. Thus the question of modular equations is 
reduced for finding a way to produce certain identities between higher dimensional 
theta functions evaluated at t and pr. In this note we apply methods from [Ko] to 
suggest such a procedure. The procedure constructs equations between 



6 



(0, r) and 9 



(0,Pt), 



for any p > 2 and yy^, are g integral characteristics. We stress that in addition to 
applications similar in the one dimensional case and possibly in cryptography we 
believe that an existence of such a procedure should lead to other applications that 
are not present in the 1-dimensional case. For example such a procedure should 
produce algebraic conditions that characterize Abelian varieties that are isogenous 
to multiplication of elliptic curves. The problem treated in this note is addressed in 
the literature with different approaches. Modular equations for hyper elliptic curves 
through I torsion subgroups of Jacobians were defined and treated in [GSj . [CKLj 
treated the case for p = 3 and [CL| seem to treat the general p along the same lines. 
Their work produces identities that use the theory of algebraic theta functions and 
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Riemann's thcta formula. The relation of these identities to the current work is not 
obvious and we plan to investigate it further in the future. 

We review the structure of this note: The first section explains the process to 
obtain identities between theta functions with integral characteristics at r and pr. 
We apply these results in second section to the one dimensional case and explain 
how this leads to a proof of existence of modular equations to the A function. This 
serves as an alternative to the classical theory of modular polynomials that has 
no analogue in the higher dimensional case. The last section we produce modular 
equations for p = 3,7. 

Acknowledgements: We thank David Lubicz whose interest in these questions 
prompted the initial motivation for this work. We thank David Kohel and Hershel 
Farkas for reading and providing valuable suggestions on an earlier version for this 
note. We especially thank the referee for very constructive remarks that substan- 
tially improved the presentation of this note. This work was partially done while 
the author visited Oakland university and the author thanks the department for 
support and hospitality. 



2. Modular equations 
We remind the reader the definition and main properties of theta functions: 

Definition 2.1. Let r be a complex g x g matrix such that: 

• T = r* i.e. T is symmetric 

• IrriT is a positive definite quadratic form. 

e 



Then t e Hg. Let 



e 



be a real 2g vector. Theta function is a complex 



analytic function on x Hg such that: 



e 



Note that the definition is the classical definition of theta function. The modern 
authors omit the factor i. We list the main properties of theta functions: 



• 6 

• 6 

• 6 



e + 2m 
e' + 2e 
e 
—e 



(z, t) — expni {e*e} Q 



(z, r) and m,e E 



(z,T) = e 



e 
e' 



(z+n+m^T^T) = exp27ri ^'"'^ — m*rrn|9 



e 
e' 



For the proof of these properties that follow from a careful series manipulation 
see [Muj , [RFj or [Koj . We remind the reader the notion of integral (rational) theta 
characteristics: 



Definition 2.2. The functions 
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e 
e' 



are called theta functions with integral (rational) characteristics if 

e Zl^ Further 



From property (2) we see that we can assume that 



e 



e 
e' 



(z, r) will be even or odd if the scalar product e'*e = 0, 1 respectively. 



This motivates the following definition: 

Definition 2.3. The integral characteristics 
0,1 



e 
e' 



is called even (odd) if e'*e 



Let us cite the duplication formula for higher dimensional theta functions that 
will be important in the sequel: 



e 

ei 



(0,r)= 



e + a 
2ei 



(0,2T)e 



a 




(0,2t) 



Here e, ei G are any g rational characteristics. The proof is [Mu| or [RF| . In 
particular replace t with 2r on both sides of the equation to obtain that: 



(4) 



e 
ei 



(0,2r)= ^ e 



e + a' 
2ei 



(0,4r)e 



a 




(0,4r) 



Let us show the following lemma: 
Lemma 2.4. Let 8 ^13^ then: 



e 



(0, r) = ^ esp {-KiE ■ /?*) e 







(0,4t) 



Proof. We write by the definition of theta functions 




e 



(0,r) = ^exp27rz|iw* + Z*^|. 

of the last equation as : 
exp2Tri | i (2m + ^0) r {2m + (if + (2m + /^)* ^ | 



Rewrite the right hand side of the last equation as 



Since exp {2T:im6) — 1 this equals to: 



exp(7ri(5*/3) x exp27ri < - 



the last sum equals to : 

Y exp {niS • e 

/3GZS 

by the definition of theta functions. 



4r m 



(3 




(0,4t) 



□ 
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Corollary 2.5. The following identity holds 



e 







(•^'4^) = ^ E exp (-^jJ* • /3) e 



S£Z!> 



(0,r) 



Proof. We can treat 



/3 




(0,4t) as unknowns in a system of linear equations 



in which the elements of the 2^ x 2^ matrix A are ea;p(7r5*/3.) Multiply A by its 
Hermitian transpose A*. Then {AA*)ii = 2^ since the diagonal element equals: 
exp(7rj(5*/3) exp(-7ri(5*/3) 29. if i 7^ j then = 0. This is because this 

element can be regarded as a sum of a nontrivial character on the group Z|. 

□ 

We apply the formula to Eq. If ei is an integer we rewrite the equation as: 



(5) 



(0,2t)= ^ exp{TTi{e + a')ei)Q 



e + a' 




(o,4T)e 







(0,4t) 



Applying Corollary 2.4 the last equation equals to : 



(6) 
and 



(o,r)e 



(0,r) 



Ce,e',a',7,5 = exp (7ri(e + Q!')*ei) X exp (7rj(e' + a')*7) x exp (7ria'*(5) 





7 



(0,r)e 



Note that for fixed (5, 7 the total coefficient of the product O 
is: 

exp (7ri(e + a')*ei) x exp (7ri(e' + a')*7) x exp {tiiol^S) 



(0,' 



which equals to unless ei + 7 + (5 = 0. In the latter case the coefficient is : 
2^ exp(7rie*5). Summarizing we obtain the following corollary: 



Corollary 2.6. 



(7) 



e 



(0, 2t) = ^ exp (7rie*<5) 6 



(o,r)e 







(0,r) 



We apply the last corollary to explain a procedure to obtain a higher dimensional 
modular equations analogues for any prime p. 



Definition 2.7. Let £)„ be a set of vectors 
• e e Z9,e, 0, 1. 



e 



such that 



2^,0 </< 2" 



For each r G Hg let a„(T) 



e 



e 



(0, t) and 



e 
ei 



ari(T) 



induces a map from 7/'„(r) : Hg ^ CP^ (I - number of vectors in Let X„ — 

ipn (Hg) . Then Xi is the image in CP^^""^ where e,ei are integral characteristics. 
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There is a map : Xn ^ Xi which omits the non integer characteristics in the 
definition of ■^n(r). 

Lemma 2.8. The map (l)n is a finite map from X„ i-^ X\. 

Proof. Because of the Transformation formula for theta functions [RF] under the 
action of Sp (g, Z) , there exists a subgroup A„ of finite index in Sp {g, Z) such that 
^„(r) induces a map /3ra(r), /?„ : Hg/A„ i-^ Xn ( Note: A„ is not a congruence 
subgroup of level n.) Hence, the map 0„ factors through a map 

:Hg/A„H^Hg/Ai, 

which is clearly finite since A„ has a finite index inside Sp {g, Z) . 

□ 

Before stating the theorem that describes the map </>„ more explicitly we state 
the following definition: 

Definition 2.9. Let if be a complex analytic domain. fi...fn : H i—> C he complex 
analytic functions. We call / constructible from /i.../n if 

• / is algebraic above C(/i.../„) 

• The Galois group of C{f) above fi-.-fn is solvable, equivalently / can be 
expressed through radical expressions involving /i.../„. 



Theorem 2.10. Let ei G Q^,eii 



2Tr. Let e e Z^. Then 6 



e 



(0,t) is con- 



structible from 6 



V 



(0, r) and rj, t]i are integral characteristics. 



Proof. Assume inductively that the theorem is true for all characteristics Q 



S 
6i 



(0,r) 



such that 5 e 1^^ and 6u 



The duplication formula implies: 



e 



(0,r)= 



e + a' 
2ei 



(0,2r)e 



a' 




(0,2r) 



But 



e 



e + a 
2ei 



(0,2t) 



satisfies the induction hypothesis. So it is constructible from 



(0,2t) 



Apply the formulas from corollary 2.6 to see that 



V 



is constructible from 



Hence 



is constructible. 



e 



e 



V 
e 



(0,2t) 
(0,r). 
(0,r) 



□ 



MODULAR EQUATIONS OF ORDER p AND THETA FUNCTIONS 



7 



Note that the proof gives shghtly more i.e. a recursive process how to construct 
the expressions 

e 



e 



from 



e 



V 



(0,r) 
(0,r) 



We reformulate it in the following corollary: 

Corollary 2.11. Let P G Xi and let Q € then there exists a process that 

produces an algebraic relationship between the coordinates of Q and coordinates of 



We rely on the last theorem and the methods developed in [Ko] to obtain process 
to produce generalized modular equation in any dimension. 

Definition 2.12. Let p be a prime. A modular equation of order p will be any 



non trivial polynomial identity between Q 
integral characteristics. 

To introduce the theorem from 
characteristics and characteristics 



(0, r) and 6 



m 



(0,pt) ?7,7yi 



Ko] we need the notion of a function order k of 
" 0...0 



0...0 

Definition 2.13. f : x Hg ^ C is an analytic function of order k and charac- 
' 0...0 ' 



teristics 



0...0 



if the following relation is satisfied: 



/ (z + n + m*r, r) — exp |27ri(— fcm*z — fcm*Tm)} / (z, t) 

Let k = P1P2 and pi,p2 are even arbitrary numbers. We quote the following 
theorem from |Ko|. 



Theorem 2.14. Let f be a function of characteristics 



0...0 
0...0 



and even order k. 



If 



1^' 



is an integral odd characteristics then the following identity is valid 



{~lY^-^' / — + T— = 



E 

and fiiy = J^i IJ^i^i- 

To obtain a modular equation for an odd number p, Choose k = 2['°3^('')]+^, I 
k — p and examine the function : 

I 



00.. .0 
00.. .0 



(z,T)e 



00.. .0 
00.. .0 



[pZ.pT) 



This function is of characteristics 



00. ..0 
00. ..0 



and order k. Define pi = 2['°92(p)l,p2 = 



2. Apply Theorem 2.14 to obtain the following: 
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Theorem 2.15. (Modular equation for p) For any 

istics: 



E (-1) 

i^,i^',0<i^'<pi,0<i^i<l 





(o,T)e 








Pi 




Pi 



odd integral character- 



(0,pt) =0 



where pi = 2^^"32{p)]_ The main theorem of the paper is given below: 
Theorem 2.16. There exist explicit equation connecting 



e 



(0,r), and Q 



(0,pr) 



XiiXi '^^6 g integral characteristics. 
Proof. According to theorem 2.15: 

Uy fi<v'^<p-L,0<Vi<l 











(o,T)e 




Pi 




pi 



(0,pt) =0 



and 



is an odd characteristics. Applying theorem 2.9 conclude O 



Pi 



(0,t) 



is constructible from O 



r] 







i 

(0,r),e 




e 


2!/. 






Pi 




Pi 



volving 6 



(0, r) and r], r/i is an integral characteristics. Replace 
(0,pr) with the corresponding radical expression in- 
(0, r) to conclude the result. 



□ 



3. The one dimensional case 



In this section we explain the connection between the theory developed in the 
last section and the usual one dimensional theory of modular equations. 

Let E be an elliptic curve given in Legendre's normal form y"^ = x{x — l)(a; — A), 
if T denotes the period that is induced by E, we have the following expression for 
A as function of r : 



(8) 

Recall the identity 

(9) 

dividing both sides by 



A(r) 



94 


" 1 " 




(0,r) 


94 


■ ■ 




(0,r) 



(0,t)=94 



(0,r) + 94 



(0,r) 



(0, r) we see that 
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(10) 

Now set e'' 
Wc write: 

(11) 
and 
(12) 



1-A(r) = 





■ ■ 
1 


(0,r) 


04 


■ ■ 




(0,r) 



(0,t) = ^o(r),e^ 



(0,r) 
OoirV 



(0,r) = ^2(r) 



The proof of the main theorem in the last section applied to the one dimensional 
case produces a homogenous radical expression of the form: 

(13) F(^,(T)^,(pr))=0. 

Divide each term of the expression by 9q{t)9q{pt). Using the definition of A(t) 
replace its quotient by A(r) and X{pt) respectively. We obtain the following classical 
theorem: 

Theorem 3.1. There exists an algebraic relation between A(t) and X{pt). 

Note that the proof of this theorem we obtain does not use the usual modular 
group theory. The proof is constructive and provides an alternative way to construct 
modular equations for any p. 

As an example consider the case p = 3 then applying the algorithm we obtain: 

(14) eo{T)eo{3T) = ^i(r)^i(3r) + ^2(r)e2(3r) 

divide the two sides of the last equation by 9o{t)9{3t) we obtain the classical 
modular equation: 



1 = ^A(r)A(3r) + ^(1 - A(r)) (1 - A(3t)). 
More details can be found in [Bo] . 

4. Examples 

We apply our theory to two cases as an example: 

4.1. p=3, g—2. In this section we outline the result of the method applied above 
to p = 3 for 5 = 2. First we observe that for g = 2 we have 6 odd and 10 even 
characteristics. We conclude immediately that overall we have 6 modular equations 
for each p. Using the Theorem 2.7 we write for p = 3 : 



Theorem 4.1. For any 

are true: 



odd integral characteristics the following identities 



J2 (_i)''^-'^''^'e 

v,v' ,0<v'<l,0<Vi<l 



(0,r)e 



(0,3r) =0 
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To achieve a more compact set of identities we rely on the classification of iden- 
tities of power 4 achieved in |AK| . 



Definition 4.2. Let 







1 



Using the last definition we can write any 2 dimensional characteristics using the 
vectors above. For example: 

' 01 



03 = 



01 



We define matrix A the encodes even characteristics in the following way : 

,(00) 



A 



11 01 10 
22 20 02 
33 21 22 



Then the following classification of 2 dimensional thcta identities is given in AK] 

Theorem 4.3. There are two types of relations for theta functions in power 4-' 

• Type I - corresponds to generalized diagonals of matrix A on the one side 



and the characteristics ( 00 
tity is true : 

^00 ~ *^11 ' *^20 

• Type II - Correspond to 2x2 sub matrices of A putting 2 not the same row 
and not the same column entries on each side of the equation. For example: 



on the other for example the following iden- 



The identities of power 4 of theta constants were obtained applying the Gauss 
elimination procedure to the matrix that defines the power 4 identities of theta 
functions 'AK] . Since the same matrix governs the identities involving terms of the 
form 0(0, r)8(0, 3r) an immediate corollary of the last theorem is the following 
classification of identities involving prime p — i: 

Theorem 4.4. There are two types of relations for theta functions involving the 
prime 3: 

• Type I - corresponds to generalized diagonals of matrix A on the one side 
and the characteristics ( 00 ) for example the following identity is true : 

eoo(0, T)eoo(0, 3r) = eii(0, r)eii(0, 3r) + 620(0, r)e2o(0, 3r) + 612(0, r)6i2(0, 3r) 

• Type II - Correspond to 2 x 2 sub matrices of A putting 2 not the same 
row and not the same column entries on each side of the equation. For 
example: 

611(0, t)6ii(0, 3r) + 620(0, r)62o(0, 3r) = 601 (0, r)6oi(0, 3r) + 622(0, r)622(0, 3r) 



This gives an algebraic way to evaluate 0y (0, 3r) from 8^ (0, t) . Compare with 
[CKL]. 
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4.2. p—7, g—2. Let us write the equations explicitly for p = 7. According to the 
recipe outlined above we need to choose k — 8 and thus our basic function will be: 



00 
00 



00 
00 



(7z, 7t) 



for p = 7 Now we use theorem 2.14 to obtain the 6 equations for each odd charac- 
teristics. For example if the odd characteristics is 



10 
10 



The equation is: 



E 

0<i^i,U2<lfi<u[,u^<2 



i-ir 



e 





(0,r)e 


1^11^2 






L 2 2 J 




L 2 2 J 



(0,7t) = 



To translate the equation into equation involving integral characteristics we use 
the duplication formulas. For example we write: 



11 

1 1 

.22. 


(0,r) = e 


■ 11 ■ 
11 


(o,2T)e 


■ 00 ■ 

00 


(o,T)+e 


■ 00 ■ 
11 


(o,2T)e 


■ 11 " 
00 



Since the other two thcta functions with integral characteristics are equal to 



0. Similar formulas hold for the other functions of the form : O 



2 2 



(0,r) 

lal to 
(0,r) 



Substituting we obtain formulas involving integral characteristics at point t, 2t, 14t. 
To reduce to equations involving r, 7r we apply the formulas from corollary 2.5 to 
functions 



e 



V\V2 
V'^V'2 



(0,2t) 



and Vi, v[ are integral characteristics. 
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